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Abstract 

We consider the application of the two-loop functional renormalization-group (fRG) approach 
to study the low-dimensional Hubbard model. This approach accounts for both, the universal 
and non-universal contributions to the RG flow. While the universal contributions were studied 
previously within the field-theoretical RG for the one-dimensional Hubbard model with linearized 
electronic dispersion and the two-dimensional Hubbard model with flat Fermi surface, the non- 
universal contributions to the flow of the vertices and susceptibilities appear to be important at 
large momenta scales. The two-loop fRG approach is also applied to the two-dimensional Hubbard 
model with a curved Fermi surface and the van Hove singularities near the Fermi level. The vertices 
and susceptibilities in the end of the flow of the two-loop approch are suppressed in comparison 
with both the one-loop approach with vertex projection and the modified one-loop approach with 
corrected vertex projection errors. The results of the two- loop approach are closer to the results of 
one-loop approach with the projected vertices, the difference of the results of one- and two-loop fRG 
approaches decreases when going away from the van Hove band filling. The quasiparticle weight 
remains finite in two dimensions at not too low temperatures above the paramagnetic ground state. 
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I. INTRODUCTION 



The discovery of high-T c superconductors, which demonstrate nontrivial properties in a 
broad temperature- and concentration range has dramatically increased interest to correlated 
low-dimensional systems, and investigation of these systems became a challenge for modern 
solid state physics. Later discovery of unconventional triplet superconductors (in particular 
Sr 2 Ru04) has attracted further attention to a possible instabilities of a Fermi-liquid state 
in the low-dimensional systems due to electronic correlations. These compounds stimulate 
theoretical interest to study the effect of correlations on electronic and magnetic properties 
of low-dimensional systems. 

The common model to treat electronic correlations is the one-band Hubbard model. At 
sufficiently large on-site Coulomb repulsion U ~ W (W is the bandwidth) this model de- 
scribes the Mott metal-insulator transition. This transition is an essentially non-perturbative 
phenomenon and is well described by the dynamical mean-field theory 1], which considers 
the limit of infinite number of dimensions and neglects spatial correlations. However, even 
in the weak- and intermediate coupling regime U < W the Hubbard model is nontrivial in 
two dimensions near some special (van Hove) band fillings or Fermi surface nesting, where 
magnetic and/or superconducting instabilities may arise [2j. The spatial correlations, not 
accounted in the dynamical mean field theory, become important in the vicinity of the cor- 
responding quantum phase transitions. Therefore, the development of methods, which are 
able to describe magnetic or superconducting fluctuations is of high interest. 

While the ID Hubbard model is exactly solvable by the Bethe ansatz and the phase dia- 
gram of this model with linearized electronic dispersion was obtained within the bosonization 
and field-theoretical renormalization group methods [3|; other numerical or approximate 
methods have to be used in higher dimensions. The applicability of numerical methods 
(exact diagonalization, quantum Monte Carlo (QMC), dynamical cluster approximations 
etc.) which treat spatial correlations is restricted by the cluster size and/or not too low 
temperatures. At the same time, there is a number of different (semi-) analytical approx- 
imations which treat the Hubbard model in the weak- and intermediate-coupling regime. 
The simplest is the mean-field approximation which treats the electron-electron interac- 

n 

tion via some effective field applied to the fermionic system [2J. Regarding the stability 
of the paramagnetic state this approach is essentially equivalent to the requirement that 
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the corresponding susceptibilities in the random phase approximation (RPA) or T-matrix 
approximation (TMA) [4| remain positive and finite in the absence of instabilities in the 
particle-hole (ph) or particle-particle (pp) channel. The corresponding electron-electron in- 
teraction vertex, irreducible in the ph or the pp channel is supposed to be equal to the bare 
on-site Coulomb repulsion U in these approaches. 

More complicated approaches account for the effect of fluctuations. These approaches 
can be subdivided into two classes: (i) approaches which consider the effect of the renormal- 
ization of the ph- or pp-irreducible electron-electron interaction vertex and (ii) approaches 
which consider in addition to (i) the renormalization of the one-particle Green functions. 
One of the approximations of the first class is the combination of RPA and TMA, which was 

n n 

proposed to account for both, the ph- and the pp-scattering [5|, [6J . In particular, one can 
use the RPA vertex (instead of the bare U) as the pp-irreducible vertex in TMA [5] , or, vice 
versa, one can use the TMA vertex, which is irreducible in particle-hole channel, instead of 
the bare vertex in RPA [6]. The two-particle self-consistent (TPSC) approximation jj] uses 
the RPA-type vertex with the effective interaction U e f instead of the bare one, the U e f is 
determined by the requirement of the fulfillment of sum rules. 

On the other hand, the commonly used approximation of the class (ii) is the FLEX 
approximation [8]. This approximaton uses the RPA interaction vertices but accounts for 



the renormalization of the one-particle Green 



of the type (ii) is the parquet approach 



unctions as well. More complex approach 
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12] which considers the contribution of 



different channels of electron scattering and their mutual interplay in the interaction vertex. 
However, the practical application of this approach for systems with the dimension ality 
d > 1 meets serious computational difficulties and was performed only in few cases [ll|, fl2| . 
The abovementioned approximations give a possibility to treat spatial correlations of the 
Hubbard model in the weak- and intermediate coupling regime. However, the accuracy of the 
results obtained within these approximations can be hardly controlled; these approximations 
can also be hardly improved. 
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The recently proposed functional renormalization-group (fRG) approaches [13| 
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19] use a different strategy. Integrating out modes with quasiparticle (qp) energy 



kk| > A, where A is the cutoff parameter, one obtains a (formally exact) hierarchy of 
RG equations for the n-particle interaction vertices. This hierarchy is usually truncated 
by neglecting higher-order vertices. To leading (one-loop) order these equations neglect 
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the six-point vertex and describe the renormalization of the two-particle electron-electron 
interaction vertices only. Therefore the one-loop fRG approach belongs to approximations 
of class (i). Unlike the RPA and TMA, however, different electron scattering channels 
are treated on the same footing within the fRG. In one dimension this approach allows to 
reproduce the results obtained earliar within the field-theoretical RG approach j2o|. The 
results for the instabili Aes^ flow of electron-electron interaction vertices, and phase diagrams 
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24j as well as the extended 25| 2D Hubbard model were 



of the standard [14i. Il6. 
also obtained at one-loop order. 

The self-energy effects, which are not included in the one-loop calculations, can be con- 
sistently taken into account at the two-loop order. In one dimension these effects are shown 

to be crucially important to describe Luttinger liquid behavior 3| . The calculation of the 
I 1 

scattering rates (l7[ , quasiparticle residues [23( and the electronic self-energy [26|, |27J in two 
dimensions using vertices obtained in the one-loop approximation showed, however, that 
contrary to the ID case the self-energy effects in 2D are much weaker. 

To estimate corrections to the one-loop approximation, however, the full calculation of 
the two-loop contributions to the flow of vertices is necessary. Contrary to the calculations 
at one-loop order, the two-loop corrections account partly for the frequency dependence of 
the vertices and their momentum dependence beyond the projection to the Fermi surface. 
Therefore, the two-loop calculations serve also as a test of the importance of the frequency- 
and momentum dependence of the vertices. Finally, they provide an information about 
quasiparticle weight, damping and the interaction-induced Fermi surface shifts. 

Although the two-loop corrections were considered previously for 2D systems in Ref. 28] 
within the field-theoretical renormalization group approach, the application of this approach 
is limited to flat Fermi surfaces and the electronic dispersion linearized near the Fermi 
surface. The advantage of the functional renormalization group approach is that it can 
be applied to both, flat and curved Fermi surfaces with or without van Hove singularities, 
since this method does not require universality of the scaling functions. The applicability of 
this approach for calculation of the two-loop corrections to scaling functions of the bosonic 



model was investigated in Ref. 



29j where the need of account of the momentum- and 



frequency dependence of the vertices was emphasized. The treatment of this dependence 
numerically is, however, a rather difficult task. 

In the present paper we use a slightly different method, which allows us to avoid consid- 
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ering momentum and frequency dependence of the higher-order vertices and calculate the 
two-loop corrections and investigate their influence on the flow of the coupling constants, 
susceptibilities and self-energies of the Hubbard model. We use the momentum-cutoff ver- 
sion of the fRG for the 1PI functions, which is applicable in the vicinity of antiferromagnetic 
or superconducting phase. 

The plan of the paper is the following. In Sect. 2 we introduce and compare the one- and 
two- loop fRG approaches. In Sect. 3 we apply the two-loop fRG approach to the ID and 
2D Hubbard models and investigate the flow of the interaction vertices and susceptibilities 
in this approach. In conclusion (Sect. 4), we discuss results of the paper and outline future 
perspectives of the method. The derivation of the two-loop equations is presented in the 
Appendix. 



II. THE MODEL AND THE TWO-LOOP FRG APPROACH 



We consider the Hubbard model 

H = - ^ tijCiaCjo + U ~ /iTl ' W 

ijcr i 

where the hopping amplitude Uj = t for nearest neighbor sites i and j and Uj = —t' for 
next-nearest neighbor sites (t, if > 0), //is the chemical potantial, corresponding to the 
particle number n. In momentum space Eq. (TjQ) reads 

H = ^£kCk CT C k(T + — ^ ^2 C kia C k 2( r' C k3<T'Ck4a^k 1 +k 2 -k3-k4, (2) 
kcr kik 2 k3k4 cr^a' 

where the Kronecker 5-symbol ensures momentum conservation, is the electronic disper- 
sion, iV is the number of sites. 

To calculate physical properties of the model ([!]) we apply the fRG approach with a 
sharp momentum cutoff (see, e.g. Ref. [23J]), which considers an effective action obtained by 
integrating out modes with the quasiparticle energy |ek| > A, A being the cutoff parameter. 
This procedure is especially convenient when there is no ferromagnetic instability developing 



in the weak-coupling regime (in two dimensions this implies if < 0.3t, see Refs. |19l. l2lj). In 
case of a ferromagnetic instability the contribution of small-momenta particle-hole scattering 
which is not included in the momentum cutoff fRG approaches may become important 



already at sufficiently large momenta scales 



m 



2l| , this case is not considered in the present 



aper. Among different versions of the 



HQ 



or one-particle irreducible (1PI) 
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HG a ppro ach (Polchinskii [l3j, Q, Wick-ordered 



23J), we use the fRG approach for the 1PI 



functions. In this approach the electron propagator at scale A has the form 

C^k,^) = — — ■ , — r (3) 



where Sa(1s, iv n ) is the self-energy at the same scale, v n are the fermionic Matsubara 
frequencies. For A > A = max(|ek|) the internal one-particle Green functions in all the dia- 
grams are zero, so that the renormalization of the physical quantities is absent: the effective 
interaction V A coincides with the bare one and T, A (k.,iu n ) = 0. The self-energy Ea(1c, iu n ) 
as well as the electron-electron interaction vertex V A (ki, fe; &3, k$) [k\, k^ and k^, foj are the 
momenta- and frequencies of the incoming and outgoing electrons, fcj = (kj, iv$)) at A < A 
can be obtained by integration of the corresponding flow equations. 

At one-loop order 5 diagrams contribute to the renormalization of the electron-electron 
interaction vertex V A {k\, k<i\ k 3 , h±) and two diagrams to the self-energy EA(k, iup) (see Fig. 
1). The corresponding flow equations can be written in the form (see Refs. [17|,li8[) 

^ = K A oS A (4a) 
^ = Va o (G? A o S A + S A o G A ) o V A (4b) 

where o is the short notation for the summation over momentum-, frequency- and spin- 
variables according to standard diagrammatic rules, see diagrams of Fig. 1, the single-scale 
propagator S\(]s.,iu n ) is defined by 

S A {k,tv n ) = — — j — - , — -. 5 

w n - e k - 0(|e k | - A)E A (k, %v n ) 

The equations (jl]) should be solved with the initial conditions V Ao = U and Ea = 0. To 
demonstrate how the fRG equations ([4]) reproduce the perturbation theory results, it is 
helpful to expand their solution in the bare interaction U. To this end, we solve them 
iteratively. Starting from the bare values of V and E we obtain after one iteration the first- 
order result for the self-energy and the second-order perturbation theory (SOPT) result for 
the vertex 

E« = UTr(G° A ) (6) 
= U + U\GloGl) 



where the index "0" stands for the bare Green functions with E = 0. After the second 
iteration we have 



E (2),l-loop = UTr ( G W ) + U 2 G A qG a qG 



A 

Ao 



vi 2)Moop = U + U 2 dA'{s£) o G« + G« o + t/ 3 [G A oG»oG»o G A ] laddcr (7) 

A 

3 /" A ° / ^00 

+U J C?A [(G A , O G A /) in O ^y(G A / O G A /) ex ] n0 n-laddor 

Here the and functions are calculated with the self-energy ([6]), "ladder" and "non- 
ladder" denote two different kinds of diagrams (see Fig. 2), "in" and "ex" denote the Green 
functions which belong to the internal and external bubble in the non-ladder diagrams, as 
shown in Fig. 2. The integrands in the second and last line do not form a total A-derivative, 
and therefore we do not obtain the exact third-order perturbation theory (TOPT) result for 
the vertex. While the integrand in the first line of the Eq. (EJ) can be changed to form the 
total derivative by the replacement Sa — > dGiy/dh, which was was shown in Ref. [3] to be 
equivalent to borrowing some terms from the two-loop corrections to the vertex, casting the 
term in the third line of Eq. ([7]) to the form of the total derivative requires full consideration 
of the two-loop corrections. 

The consideration above provides a definition of the n-loop approximation as an approx- 
imation which correctly reproduces n-loop parts of the diagrams for the two-particle inter- 
action vertex and n + 1-loop parts of the self- energy diagrams. In the presence of logarithic 
divergencies (e.g. in one dimension), when G A oG A ~ ln(A/a) (a <C A is a small parameter), 
the terms {G% o G% o G%) o (dG% / dA') and U 3 {G° A , oG A ,) cx od{G A , oG A ,) in /c/A' which appear 
at the two-loop order (see below) and are necessary to combine to the total derivatives in 
the Eq. ([7]) can be neglected to leading logarithmic order. This provides another, more 
conventional definition of the n-loop approximation as an approximation which correctly 
treats the terms U m \n m ~ n (A/a) in the perturbation series for the vertex (m > n). Note, 
however, that in two dimensions the bubbles G A , o G A , are either non-divergent for arbitrary 
fillings or may contain squared logarithmical divergencies for some special (van Hove) band 
fillings, and the latter definition of the n-loop approximation can not be applied. 

To go beyond the one-loop order of the Eqs. (j4j) one has to take into account the 
contribution of the three-particle interaction vertex (see Refs. 17J, llSj and Appendix). 



Generally, this vertex generates contributions to the two-particle interaction vertex V with 
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an arbitrary number of loops n > 3. At this stage two different approximations are possible, 
(a) keeping only contributions which are necessary to treat exactly diagrams with fixed 
number n of loops and (b) keep all the contributions which are generated by an integration 
of the equation for the n + 1-particle vertex, neglecting n + 2-particle vertex. In the present 
paper we restrict ourselves to the approximation (a), i.e. consider only those contributions 
to RG flow which are necessary to treat exactly the two-loop parts of the diagrams. 

At the two-loop order the flow equation for the self-energy (j4aj) does not change while 
32 new diagrams contribute to the renormalization of the vertex (see Fig. 3). The resulting 
two-loop equations have the form (see Appendix for the derivation) 



dA 
dV A 



V A o S A (8a) 
V A o (G a o S a + S A o G A ) o V A (8b) 



dA 

Ao 



+S A o J dA'V A / o G A i o V A i o G A / o V A i o S A > 



A 



Contrary to the one-loop approximation the frequency dependence of the vertex becomes 
essential at the two-loop order. This can be seen from the fact that to reproduce the 
TOPT result one needs to iterate Eq. fl8bl) twice. If one neglects the frequency dependence 
of the vertices, the equations ([8]) fail to reproduce correct TOPT results. The necessity 
of the account of the frequency and momentum dependence of the vertex was previously 
emphasized in the two-loop calculation of t 
calculation in the 2D Hubbard model 



22 



re /9-function of <p theory [29] and the self-energy 



23 



26. 



27). 



To avoid having explicitly the frequency- and momenta- dependent vertices, we integrate 
Eq. (18bl) formally and keep frequency- and momentum dependence coming from the one-loop 
term only to obtain 



Va 


= v A + SV A 




Ao 


6V A 


= J dA'[V A , 




A 



(9) 



where V a = VVa and the operator V projects the external frequencies to zero and the exter- 
nal momenta to the Fermi surface. Reinserting this vertex into the one-loop contributions 
of Eqs. (El) and using projected vertices in the two-loop contributions, we obtain to linear 
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order in 5V 

Ao 



— — - = Va o Sa + Sa ° / dA'[V A > ° Ga< ° Sa> ° V A * — V(V a> ° G A i ° S A < ° Ka/)] (10a) 

A 

= V [V A o (G A o S A ) oV A -V A o (G A o S A ) o (Vf -V ) 



Ao 

+Va ° (Ca ° S'a) ° J dA'V A > o Ga' o S'a' o V A > 



A 

Ao 



+S A o y dA'V A > o G A / o V A , o Ga 1 o V a > o S a > | (10b) 



where 



Ao 

Vf = V + V J dA'(VA' o G A i o S A ' o V a > 



A 



is the analogue of the 1-loop vertex calculated with the two- loop vertices V a>- After this 
reinsertion, only the projected vertices V enter the Eqs. ( TlOl) . While the last term in the Eq. 
( 110b!) accounts for the two-loop corrections to the flow, other integral contributions correct 
the effect of the vertex projection in one-loop diagrams. In particular, for A-independent V 
and £ the last two terms in the Eq. fllObl) combine to a A-derivative of the corresponding 
two-loop diagram, so that these equations with the initial condition v\ o = Va = V = U, 
S Ao = correctly reproduce the result of the TOPT after one iteration. The two-loop 
fRG equation for the self-energy fllOal) is identical to that investigated earliar with one- loop 
vertices j^jj • The flow of the susceptibilities is described by the equation, similar to the Eq. 

(HQED, 



dx A 



T A o (G A o S A ) o T A (11a) 
-A = v[T A o(G A oS A )oV A -T A o(G A oS A )o(V A L -V ) 



Ao 

+T A o (Ga ° Sa) ° J dA'V A i o G A > ° S'a' ° V a> 



A 

Ao 



+S A o J dA'T A ' o G A i o Vm ° G\> ° V A > o S A > | (lib) 

with the initial condition x Ao = an d T Ao is determined by the symmetry of the order 
parameter, e.g. Ta = 1 for the antiferromagnetic and singlet superconducting susceptibility, 



T\ (k) = cos k y — cosk x for the d-wave superconducting susceptibility etc. 

Let consider the local in A version of the equations ffTUl) . which is obtained by the replace- 
ment V\i — > V\, Ea' — > S A . This replacement introduces corrections of the same order, 
which are already neglected in the two-loop approximation, and therefore, can be considered 
on the same level of an approximation. In this way we obtain 



dA 
dV A 
dA 



V A o S a + S A o V a 



[1-V){GaoGi 



oVa 



(12a) 



V{V K o A( Ga oGa)oFa 



+V a o (G a o S a ) o V a o [(1 - V)(G A o G A ) 
+5a ol/ A oG A ol/ A oG A oy A o Ga} 



oV^a 



(12b) 



The local equations (|T2|) have a similar form as the two-loop equations in the field-theoretical 
approaches, e.g. for the ID fermionic systems [3j], and, therefore, can be used to make 
connection with these approaches. The terms with the projection operator V coming from 
Eq. substract the one-loop (In 2 in ID case) contributions from the third-order diagrams 
for the vertex. The coresponding contribution to the self-energy (last term in the first 
equation) is a frequency- independent constant, which can be omitted. Contrary to field- 
theoretical approaches, the Eqs. f|T2l account for both, regular and singular terms in the 
perturbation expansion and are written for the coupling constants themselves, not for their 
invariant combinations with the self-energy. 

To solve numerically Eqs. (11 01) . ffTTl) . or (1121) . we use the discretization of momentum 
space in 2 patches (L and R) in one dimension and N p = 24 patches in two dimensions with 
the same patching scheme as proposed in Ref. 19]. For the frequency dependence of the 
self-energy we use N u = 100 frequencies Wi suitably chosen on the imaginary axis (these 
frequencies do not have to coincide with the Matsubara frequencies since for a frequency- 
independent V the self-energy is defined on the entire imaginary axis, cf. Ref. |26|). We 
account for the self-energy effects by expanding the self-energy Sa(w) around v = and 
introducing Z-factors 



[i - dim£ A (Wv)/cH;;i 



o- 



(13) 



the feedback of the imaginary part of the self-energy to the flow of vertices is neglected. The 
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Green functions in this approximation take the form 23 j 

kk|-A) 



G A (k,iv n ) 
S A (k,iu n ) 



'A 



(14) 



ZtS(\e k \-A) 



w 



n 



where kp corresponds to the projection of the vector k to the Fermi surface. 

The derivative of the self-energy which enters the Eq. (|T3l) is determined numerically from 
the values of the self-energy at the first two frequencies iuip- The approximation ( 1T41) can be 
applied only in the paramagnetic state without strong exchange and/or umklapp scattering 
(i.e. away from half filling in one dimension and at not too low temperatures and not too close 
to the van Hove band fillings in two dimensions). Above the antiferromagnetically ordered 
ground state the divergence of the corresponding vertices leads to a pseudogap structure 



of the self-energy and spectral functions 



26, 



271 ] . This structure can be correctly described 



only with the frequency-dependent self energy and is not considered here. We also neglect 
the first and third terms in the flow equations for the self-energy ( fTUj) and ( TT2T) as responsible 
purely for the deformation of the Fermi surface by the interaction. This deformation was 
found numerically to be small in two dimensions at small next-nearest 



can be treated accurately by introducing corresponding counterterms 



lopping t' [17J and 



15 



22, 



331 . 



III. RESULTS 



A. ID case 



First we consider the results for the ID electronic dispersion. 



Ek = — 2t cos k — fi 



(15) 



In this case we have only 2 patches (L and R) at kp = ± arccos(/x/2). After the pro- 
jection to the Fermi points, only 4 independent vertices remain: V\ = V(L, R; R, L), 
V 2 = V(L,R;L,R),V 3 = V(L,L;R,R), and V 4 = V(L,L;L,L) = V(R,R;R,R). With 
the linearization of the dispersion (|15p near the Fermi points, the flow of these vertices 
in the two-loop approximation is well-studied in the field-theoretical approach [3| and it is 
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described by the equations: 



dg 3 /dl 



dg^/dl 



dg' 2 /dl 



dgi/dl 




y j 1 a 9 2 
(929l ~ 9 3 i) 



9AW2) 2 + ai - 29*94] 



(16) 



Aix 2 v F 



d\rt Z/dl 




t 2 1 2 1 2\ 

(#1 - 9192 + 92 +93) 



(17) 



where g>j = Z 2 Vi are the invariant coupling constants, g' 2 = g± — 2g 2 , vp = 2t sin hp, Z = In A. 
We emphasize once more, that the difference of the equations ( f!2l) and ( flBT) is that the latter 
account for the universal contributions to the flow of the coupling constants only, while the 
former treat also the non-universal contributions, e.g. connected with the nonlinearity of 
the dispersion. 

The result of the solution of Eqs. (fl2|) for gi = t, g 2 = 2t, g 3 = g± = \i = 0, and T = lCT 4 t 
is presented in Fig. 4. Chosing g 3 = allows us to avoid the effect of umklapp scattering. 
To verify that the result of the solution of the Eqs. ( Till is indeed reproduced at A C 1, 
we use the result of integration of Eqs. (112]) at some cutoff parameter A c = e~H < 1 as 
a starting condition for the Eqs. (TIB"]) and compare the result of the solution of the Eqs. 
(fl2"|) and (TTBT) for A < A c . One can see that the results of the local two-loop fRG approach 
agree with the corresponding results of the field-theoretical approach, Eqs. ffT6|) . At the 
same time, the results of the solution of local and nonlocal fRG equations are different due 
to the nonuniversal initial part of the flow. We have verified that this is mainly connected 
with the momentum dependence of the Fermi velocity, the difference almost disappears for 
the linearized version of the dispersion (fl5l) . The nonlocal equations fllOap are expected to 
treat better the effect of the nonlinearity of the dispersion, therefore we consider only their 
solution in the 2D case below. 

B. 2D case 

For the discussion of the results of fRG approach in 2D case we also consider the solution 
of the Eqs. ( flOi) without the two-loop corrections (i.e. without the last term in the second 
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equation) to investigate how the one-loop flow changes due to correction of the errors of 
vertex projections by the second and third term in the Eq. fllObl) . The difference of the 
latter results from the one-loop results shows the effect of the vertex projection on the 
renormalization group flows, while their difference to the two-loop results shows the effect 
of the two-loop corrections. 

First we consider the results of the solution of Eqs. (flOj) for the dispersion 

£k = — t(cosp x + cospy) + t\ cosp x — cospy\ — /i (18) 



28 



34j. 



The corresponding Fermi surface has flat parts along the directions \p x \, \p y \ = arccos 
The field-theoretical approach for a flat Fermi surface was applied earliar in Ref. 
The results of the numerical solution of Eqs. fflUj) for fi — t, U — 7.81t are shown in Fig. 
5. We choose this relatively large value of the interaction U since it corresponds to the 



value of the dimensionless coupling constant UA/(tcvf) = 3 used in Refs. |28|, |34j where 
vf is the Fermi velocity, A is the length of the Fermi surface flat parts (vf = a/4£ 2 — /i 2 , 
A = 2 arccos(— /// (2t)) for the dispersion ({TBI ); a larger value UA/ ijivp) = 10 was considered 



in Ref. 



28] . One can see that the vertices without the two- loop corrections diverge at much 



larger energy scales compared to the two-loop results in agreeement with Refs. [28l. |34J|. The 
scale of the vertex divergence in the one-loop calculation with partly corrected projection 
errors (Eqs. (flOl) withot the last term in the second equation) agrees with the result of the 
one-loop field-theoretical approach (not having any projection errors), but is larger than the 
scale of the vertex divergence in one-loop approach with vertex projection. The flow of the 



vertices in the two-loop approach agrees also with the results of Refs. [28l. 134 ] . 

The functional renormalization group approach can be further applied to the t-t' Hubbard 
model with the dispersion 

£k = — 2t(cosp x + cospy) + 4t'(cosp x cospy + 1) — ji (19) 

where the conventional field-theoretical approach is not applicable due to the presence of 
squared logarithmic singularities in the perturbation series near van Hove band filling (fi = 
0). The results for the flow of vertices and susceptibilities for t' = O.lt, U = 2t and the 
fillings close to vH band filling (fi = ±0.1t) are presented in Fig.6. At p = O.lt (above vH 
filling) the largest susceptibility is observed with respect to the antiferromagnetic instability 
in both, the one- and the two-loop approaches. For this value of fi and chosen temperature 
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T = O.lt the results for the susceptibilities in one- and two-loop approach substantially differ. 
The antiferromagnetic susceptibility in the end of the flow in the one-loop approach with 
partly corrected errors of vertex projection (Eqs. (Ill) I) with omitted last term in the second 
equation) is larger than the results of this approach with vertex projection and deviates more 
from the two-loop results. Therefore, the results of one-loop approach with vertex projection 
agree better with the two-loop results, which is possibly due to account of only universal 
terms of the flow in these approaches. For fi = — 0. It (below vH filling) we observe the same 
qualitative tendencies with smaller difference of the results of one- and two-loop approaches. 
With decreasing temperature, the superconducting instability becomes dominating in this 
case (see below). 

The results for the flow of the vertices and susceptibilities at the filling further from vH 
one (/i = — 0.5i) are shown in Fig. 7. At this filling and not too low temperatures the anti- 
ferromagnetic susceptibility is the largest one (not shown), but with decreasing temperature 
the d-wave superconducting instability becomes the leading instability. Susceptibilities in 
the one-loop approach with partial correction of vertex projection errors are larger than the 
one- and the two-loop approaches. The susceptibilities in the end of the flow of one- and 
two-loop approaches are, however, closer to each other than in the above considered case 
li = -0.lt. 

The calculated temperature dependences of the susceptibilities for antiferromagnetic and 
superconducting instabilities, as well, as Z-factors for fi = ±0.1t are shown in Fig. 8. 
At n = O.lt we observe a maximum of the antiferromagnetic susceptibility in the two- 
loop approach, while the corresponding susceptibility in the one-loop approach diverges 
with decreasing temperature (Fig. 8a). More generally, we find that the divergence of 
the vertices (and susceptibilities) is strongly suppressed in the two-loop approach. This 
divergence is not, however, fully removed, since for smaller /i < 0.08 we find again the 
possibility of the antiferromagnetic ground state in the two-loop approach. At fi = —O.lt 
the antiferromagnetic susceptibility also has a maximum at some temperature and then 
decreases with decreasing T (Fig. 8b), while the superconducting susceptibility increases, 
showing the possibility of the superconducting ground state in both, one- and two-loop 
approaches. The increase of Xdsc m the one-loop approach is again more pronounced than 
in the two-loop approach, so that the temperature where the susceptibility diverges in the 
two-loop approach is expected to be much smaller than in the one-loop approach. The 
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Z-factors decrease almost linearly with In (t/T) at intermediate temperatures, but below 
the temperature, where the maximum of the susceptibility is reached, their temperature 
dependence becomes linear in T and therefore Z kp are not expected therefore to vanish at 
lower temperatures. We have also verified during the calculations that the imaginary part 
of the self-energy remains linear function of v at small imaginary frequencies. 

In Fig. 9 we summarize the results for the qp damping Tk F = — lm£(kp,0) and the 
FS shift Re£(kp,0) estimated at different FS points in the end of the two-loop fRG flow. 
The qp damping depends almost linearly on temperature at not too low temperatures; this 
dependence becomes quadratic at low T. The observed linear dependence of the scattering 
rates at not too low temperatures may be due to closeness to the antiferromagnetic quantum 
critical point; more detail investigations of this dependence are, however, required. At the 
same time, the quadratic temperature dependence at low temperatures supports the Fermi- 
liquid picture in this temperature range above the paramagnetic ground state. The Fermi 
surface shift contributions are small and negative for [i < 0. For /i > the Fermi surface 
shifts have opposite signs at the point closest to the (ir, 0) point and to the diagonal, leading 
therefore to small deformation of the Fermi surface, which makes it flatter. 

IV. DISCUSSION AND CONCLUSIONS 

We have considered the effect of the two-loop corrections on the fRG flow. In ID case we 
find that the nonuniversal corrections contribute to the flow at large momenta scales, while 
at small momentum scales we have recovered the results obtained from the field-theoretical 
RG approach. For the 2D case with flat Fermi surface we also find good agreement with 
the previous results of the field-theoretical RG. The fRG approach was applied further to 
the case of curved Fermi surface without nesting, where we obtained the flow of the vertices 
and susceptibilities at the two-loop level. 

In two-dimensions the two-loop corrections do not change the leading instability, but may 
lead to slight shift of the phase boundaries in comparison with the previous one-loop analysis. 
The difference of the two-loop results and the one-loop results with projected vertices in two 
dimensions is smaller than to the one-loop results with partly corrected projection errors 
and decreases going away from the van Hove band filling. Therefore, the commonly used 
one-loop fRG approach with projected vertices serves as a good starting point for calculating 
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higher-loop corrections. 

We have also considered the flow of the qp spectral weight, qp damping and the Fermi 
surface shift. In agreement with earliar studies, in the considered cases of not too low 
temperatures in two dimensions we obtain the qp weight Z ~ 0.9, so that the quasiparticles 
remain well defined during the fRG flow. The qp damping and estimated Fermi surface 
shifts are also numerically small. 

Possible future applications of the method would be its implementation within the 
temperature-cutoff fRG scheme [3], where the two-loop corrections are expected to be 
smaller than for momentum cutoff due to better treatment of degrees of freedoom with 
different excitation energy. The calculation of the two-loop corrections for the temperature 
cutoff is, however, more difficult task since it requries more intense numerical calculations 
caused by the smoothness of the cutoff. Another possible extension of the method would 
include consideration of the frequency dependence of the self-energy and/or vertices, which 
also has to be performed. 
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VI. APPENDIX. DERIVATION OF THE TWO-LOOP FRG EQUATIONS 

In this Appendix we consider the derivation of the two-loop RG equations. We use the 
notations of Refs. [pr], Q, which considered the 1PI RG equations for the terms of the 
expansion of the 1PI generating functional in fermionic field, 

m>0 

where <fi{X) are fermionic fields, 

X = (x,r,a,±) (20) 

is the short notation for the space, time, spin, and charge variables (the ± sign corresponds to 
the incoming and the outgoing particles, respectively). At the two-loop order the hierarchy 
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of RG equations for the 1PI functions is truncated at the three-particle vertex and has the 
form Q Q 

7 i 2) (0) = (0,Q0) + ^Tr[S A 7i 2) ] 

7i 4) (0) = ^Tr[S A 7i 4) ] - ^Tr[5 A 7i 2) G A 7f] 

7a%) = -^Tr[S A 7i 4) G A 7i 2) + S A ^G A ^} + ^[Sa^Ga^Ga^] + 0( 7 f)(21) 
where 

the trace is taken with respect to X-variables and dots denote derivatives with respect to 
A. 

For practical calculations the vertices 7^ (0) are expressed explicitly through the 0-fields 

as 

7i m) (0) = / ^ m 7 m (A|X)0(X 1 )...0(X m ) (23) 

ml z — ' J 

m 

and 

7i m) (X,y,0) = — I d m X' lm+2 (A\X,Y,X')<p(X[)...<p(X'J (24) 

771. j 

In these notations the RG equations ( 12T1) read 

72 (A|X) = i y d 2 Y l4 (A\X,Y)S A (Y) (25a) 
74 (A|X) = 1 J d 2 Y l6 (A\X,Y)S A (Y) - l - J d 4 YB A (X, Y)L A {Y) (25b) 
7 6 (A|X) = iy d 6 rD A (X,F)M A (F)-i y d 4 y£ A (X,r)L A (T) (25c) 



where 



£aOO = Sa(Y 1s Y 2 )G A (F 3 , Y 4 ) + G A (Y U Y 2 )S A (Y 3 , Y 4 ) (26) 
5 A (X,r) = 1a (A\X 1 X 2 ;Y 2 ,Y 3 ) 14 (A\X 3 X 4 ;Y,,Y 1 ) 
^(AlXxXs; y 2 , r 3 )74(A|X 2 X 4 ; Y 4 , Yi) 

+ 74 (A|X 1 X 4 ; F 2 , y 3 )74(A|X 2 X 3 ; F 4 , Yi) (27) 
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and 

M A (Y) = S A (Y h Y 2 )G A (Y 3 ,Y 4 )G A (Y 5 ,Y 6 ) 
+G A (F 1 , Y 2 )S A (Y 3 , Y 4 )G A (Y 5 , Y 6 ) 

+G A (Y h Y 2 )G A (Y 3 , Y 4 )S A (Y 5 , Y 6 ) (28) 
D A (X,Y) = 14 (A\X 1 ,X 2 ;Y 2 Y 3 ) 1a (A\X 3 X 4 ;Y a Y 5 ) 1a (A\X 5 X 6 ;Y 6 Y 1 ) 

+14 permutations (29) 
E A (X,Y) = 7 4 (A|X 1 X 2 ;r 2 F3)76(A|X 3 X 4 X 5 X 6 ;r4^i) 

+14 permutations (30) 

As discussed in main text, in the present paper we consider the expansion in the number 
of loops rather than in effective vertices. Iterating the Eq. (I25cl) for 7W, one can easily 
see that the first term in this equation corresponds to the three- loop contribution to 7*- 4 - 1 . 
Therefore, in the following we neglect this term. Due to this neglection, the Eq. ( 125c[) can 
be integrated analytically. Substituting the result into Eq. (I25bl) we obtain 

72 (A|X) = \j d 2 Y l4 {A\X,Y)S A {Y) 

Ao 

<y 4 (A\X) = \J dM J d 2 Y J d 6 Y'S A (Y)M A (Y')D A ((X,Y),Y') 

A 

~J d 4 YL A (Y)B A (X, Y) (31) 

This substitution considerebly simplifies the numerical solution of the equations since one 
has to consider the vertices 7 2 and 7 4 only. 

Introducing the self-energy and 2-particle irreducible vertex by 

£(£i,£ 2 ) = 7 2 ((£i,T,+),(£ 2 ,T,-)) 
V(^ 2 ; U*) = 7 4 ((£i, T, +), (£2, I, +), « 8 . T> -)> (£4, I, -)) (32) 

where £ f = (xj, n) and the Fourier transformed quantities 

E(k) = J d^^Qe^-^ 
V(hk 2 ;k 3 k 4 ) = Jd^V(^2-,Ud^ lkl+ii2k2 ~ ii3k3 ^ ik4 (33) 
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where fcj = (k,, iun ) arid exploiting charge-, spin- and translational invariance in the same 



way as in Refs. 
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18| we obtain the equations (1251) of the paper. 



181 ] we introduce 



Similar derivation can be performed for susceptibilities. Following Ref. 
vertices <-y^ n,?1 * ) where m refers to the number of boson lines and n to the number of fermion 
lines, which enter or go out of the vertex. The equations for the vertices ^^ n,ri ^ read 



7i M) (0) 



(0,2)i 



(34) 
(35) 

(36) 



Expanding again •y( m ' n \(p) in 0, neglecting first three terms in the equation for 7a' 4 '' (0) as 
corresponding to the higher-loop order, substituting the result for 7 A 1 ' 4 ' ) (0) in the equation 
for 7 A 1,2 ' ) (0), we obtain 



7 (1 ' 0) (A|X, X') 
7 (1 ' 2) (A|X,X') 



rf 4 FL A (F) 7 (1 ' 2) (A|X, Y 2 , y 3 )7 (1,2) (A|X', F 4 , Y x ) 

Ao 

dM ! d 2 Y I d e Y'S A (Y)M A (Y')D A (X,(X',Y)X) 



(37) 



-\ I <l l YL A (Y)B A (X.X'.Y) 



(38) 



where 



B A (X,X',Y) 
D A (X,X',Y) 



7(i, 2 )(A|X; Y 2 , Y 3 ) l4 (A\X[X' 2 ; F 4 , Y 1 ) (39) 
7(i, 2 )(A|X; Y 2 Y 3 ) l4 (A\X[X' 2 ; Y,Y 5 )^(A\X^X' 6 Y 6 Y) 
+permutations (X r ) (40) 
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FIG. 1: The diagrams for the self-energy (a) and vertex (b) flow at the one-loop order. The solid 
lines correspond to the cut propagator G\, the lines with dash - to the single-scale propagator S\, 
boxes - to the vertices V\. Lines inside the box show the direction of spin conservation. 



FIG. 2: Ladder-type (left) and non-ladder (right) diagrams in the third order perturbation theory. 
Solid lines correspond to the bare electronic prpagator Go, boxes - to the bare interaction U, "in" 
and "ex" denote internal and external bubble in the non-ladder diagrams. 

[32] A. Neumayr and W. Metzner Phys. Rev. B 67, 035112 (2003) 
[33] M. Salmhofer, Ann. Phys. 16, 171 (2007). 
[34] H. Freire, private communication 
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FIG. 3: The diagrams for the contributions to the flow of the vertex at the two-loop order in the 
particle-particle (a) and the particle-hole (b) channels. The two of the three lines without dash 
correspond to the G\i propagator and one to the S\i propagator (the circle arrow with double 
dash denotes their permutations), the line with dash corresponds to the single-scale propagator 
S\. The other notations are the same as in Fig. 1. 
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FIG. 4: (color online) The flow of the vertices (a), the antiferromagnetic (AF) and the singlet 
superconducting (SSC) susceptibilities (b) of the ID Hubbard model within the non-local (solid 
lines) and local (dot-dashed lines) two- loop approaches and the one-loop approach (dahed lines). 
The results of the solution of the field-theoretical two- loop equations ( [T6|) with initial vertices, 
obtained in the fRG approach at the scale A = e~ A t (marked by arrow) are shown by bold lines 
with circles. 



23 





1 




• i 




80 


1L 
q1L 

2L 

— 2Lz 


)2 ! g2 j 
I 

i ! 




- 


60 




ii 






40 


- 


i I 91 !' 

I' 
II 

) I !' 

! f 
> 




■ 


20 


/' 

/' 

/ 

./ i 


f 

/.'/ | 

/// // 

-y^y 




gi 


n 


i 




i 





2 4 6 

ln(1/A) 



FIG. 5: (color online) The flow of the vertices gi = 7(1,12,12) and g 2 = 7(1,12,12) of the 
2D Hubbard model with the flat Fermi surface in the 24-patch one-loop approach with vertex 
projection (1L, dashed lines), the one-loop approach with partly corrected errors of the vertex 
projection (qlL, dot-dashed lines) and the two-loop approach with account of Z-factors (2Lz, solid 
lines) and without Z-factors (2L, dash-dot-dot lines). The first and twelth patches being close to 
the centers of the opposite FS sides. 
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FIG. 6: (color online) The flow of the vertices g x = V(l,7,7), g 2 = V(l,7,l), g 3 = V(l,l,7), 
and 34 = V(1 3 1,1) (a,c), the antiferromagnetic (AF), d-wave superconducting (dSC) and the 
ferromagnetic (F) susceptibilities (b,d) of the 2D t-t' Hubbard model with U = 2t, t'/t = 0.1, 
T = O.li and fi = O.lt (a,b), \x = — 0.1t(c,d) in the 24-patch one- and two-loop fRG approaches. 
The first and seventh FS patch correspond to points, closest to two different van Hove singularities. 
The solid lines in b) and d) correspond to the two-loop approach; the dashed lines - to the one- 
loop approach with projected vertices, the dot-dashed lines - to the one-loop approach with partly 
corrected errors of the vertex projections, other notations are the same as in Fig. 5. 
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FIG. 7: (color online) The flow of the vertices (a) and the susceptibilities (b) of the 2D t-t' 
Hubbard model, U = 2t, t'/t = 0.1, \x = — 0.5i, and T = 0.025t in the 24-patch one- and two- loop 
fRG approaches, the notations are the same as in Fig. 6. 
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FIG. 8: (color online) The temperature dependence of the antiferromagnetic and d-wave supercon- 
ducting susceptibilities in one- (dashed lines) and two-loop (solid lines) functional renormalization 
group approaches to the 2D t-t' Hubbard model with U = 2t, t'/t = 0.1, jj, = O.li (a) and n = — O.lt 
(b). The temperature dependence of the Z-factors (right axis) is shown by dot-dashed lines; the 
number of dots corresponds to patch number. 
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FIG. 9: (color online) The temperature dependence of the qp damping Tk F = — ImE(kp,0), (a,b) 
and the FS shift ReX(kp,0)/t (c,d) at different FS patches (1,2,3) obtained in the 24-patch two- 
loop fRG approach at U = 2t and t'/t = 0.1, /j. = O.lt (a,c) and /i = -0.lt (b,d). The patch 1 is 
the closest to the (ir, 0) point. 
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